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ABSTRACT 

We present new optical and Chandra observations of the field containing the ultraluminous X-ray 
source NGC 1313 X-2. On an ESO 3.6 m image, the Chandra error box embraces a R = 21.6 point-like 
object and excludes a previously proposed optical counterpart. The resulting X-ray/optical flux ratio of 
NGC 1313 X-2 is ~ 500. The value of fx/ f opt, the X-ray variability history and the spectral distribution 
derived from a re-analysis of the ROSAT, ASCA and XMM data indicate a luminous X-ray binary in 
NGC 1313 as a likely explanation for NGC 1313 X-2. If the X-ray soft component observed in the XMM 
EPIC spectrum originates from an accretion disk, the inferred mass of the compact remnant is ~ 100M Q , 
making it an intermediate mass black hole. The derived optical luminosity (L w 10 5 Lq) is consistent 
with that of a w 15 — 20M Q companion. The properties of the environment of NGC 1313 X-2 are briefly 
discussed. 

Subject headings: galaxies: individual (NGC 1313) — stars: individual (NGC1313 X-2/MS 
0317.7-6647) — X-rays: binaries — X-rays: galaxies 



1. INTRODUCTION 

First revealed by the Einstein Observatory (see e.g. 
Fabbiano 1989), point-like, off- nuclear X-ray sources with 
luminosities significantly exceeding the Eddington limit for 
one solar mass are being progressively discovered in the 
field of many nearby galaxies. To date, hundreds of such 
sources have been found in dozens of galaxies, both el- 
lipticals and spirals (e.g. Colbert & Ptak 2002). These 
powerful objects, commonly referred to as ultraluminous 
X-ray sources (ULXs), do not appear to have an obvious 
Galactic counterpart. Despite some of them have been 
identified with supernovae or background active galactic 
nuclei, the nature of most of these sources remains un- 
clear. X-ray spectra have been obtained with the Einstein, 
ROSAT, ASCA and recently XMM-Newton and Chandra 
satellites. Although the statistics is rather poor, in many 
cases fitting with simple models indicates that the spectral 
properties are consistent with those of Galactic black hole 
binaries (e.g. Foschini et al. 2002a). About half of them 
show some degree of variability in the X-ray flux. Among 
the various possibilities, the most favored explanation is 
that ULXs are powered by accretion and that they are 
somewhat special X-ray binaries, either containing an in- 
termediate mass black hole (BH) with Mbh ^ 100 Mq 
(e.g. Colbert & Mushotzky 1999; Kaaret et al. 2001) or 



having beamed emission toward us (e.g. King et al. 2001; 
Kaaret et al. 2003). For a recent review on the properties 
of ULXs we refer to Fabbiano & White (2003). 

Optical observations are of fundamental importance to 
better assess the nature of these sources but they are still 
rather scarce (see e.g. Cagnoni et al. 2003; Foschini et al. 
2002b) . Some ULXs have optical counterparts in the Dig- 
itized Sky Survey or Hubble Space Telescope images (e.g. 
NGC 5204 X-l; Roberts et al. 2001; Goad et al. 2002) and 
some appear to be embedded in emission nebulae a few 
hundred parsecs in diameter (Pakull & Mirioni 2002). 

NGC 1313 X-2 was one of the first sources of this type 
to be found. It was serendipitously discovered in an Ein- 
stein IPC pointing toward the nearby SBc galaxy NGC 
1313 (Fabbiano & Trinchieri 1987). Originally included 
in the Einstein Extended Medium Sensitivity Survey as 
MS 0317.7-6647, it is located ~ 6' south of the nucleus of 
NGC 1313. Stocke et al. (1995) investigated the nature 
of MS 0317.7-6647 on the basis of X-ray, optical and ra- 
dio observations. They identified a possible optical coun- 
terpart and concluded that the source could be either a 
Galactic isolated neutron star or a binary containing a 
massive BH in NGC 1313. Spectral fits to ROSAT PSPC 
data (Stocke et al. 1995; Colbert et al. 1995; Miller et al. 
1998) yielded results consistent with many single compo- 
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nent models. ASCA observations (Petre et al. 1994; Mak- 
ishima et al. 2000) are described successfully by a multi- 
color disk blackbody (MCD) model, representing thermal 
emission from a standard accretion disk around a BH. A 
very recent analysis of a XMM EPIC-MOS observation of 
NGC 1313 (Miller et al. 2003) indicates that two spectral 
components, soft and hard, are required to fit the spec- 
trum of NGC 1313 X-2 and the normalization of the soft 
component yields a conspicuous mass of the black hole 
M BH > 830M Q . 

We present new optical 2 and Chandra observations of 
NGC 1313 X-2, with the aim to shed further light on its 
enigmatic nature. In § 2 we present Chandra data and re- 
analyze all the available X-ray observations of NGC 1313 
X-2. In § 3 optical observations of the field of this ULX 
are reported. Finally, the implications of our results on 
the nature of NGC 1313 X-2 are discussed in § 4. 

2. X-RAY DATA 

NCG 1313 X-2 was first observed by Einstein with the 
IPC instrument in 1980. It was then pointed several 
times by ROSAT (PSPC and HRI) between 1991 and 1998 
(Stocke et al. 1995; Colbert et al. 1995; Miller et al. 1998; 
Schlegel et al. 2000), by ASCA (SIS and GIS) in 1993 
and 1995 (Petre et al. 1994; Makishima et al. 2000) and 
by XMM (EPIC) in 2000 (Miller et al. 2003). Most recent 
data are from a 2002 Chandra (ACIS-S) pointing, reported 
here for the first time. The Chandra observation began on 
13 Oct 2002 and had a duration of 19.9 ks. The primary 
goal of the observation was to study sources near the cen- 
ter of the galaxy, but the aim-point was adjusted to also 
place NGC 1313 X-l, NGC 1313 X-2, and SN 1978K on 
the S3 chip of the ACIS-S. In addition, we present here a 
complete analysis of the XMM observation, including the 
EPIC-PN data, which were not considered by Miller et al. 
(2003). 

2.1. X-ray astrometry 

An accurate determination of the X-ray position of NCG 
1313 X-2 can be obtained from the 2002 Chandra pointing, 
using the Chandra aspect solution. Chandra data were ex- 
tracted from the S3 chip on the ACIS-S and subjected to 
standard processing and event screening. No strong back- 
ground flares were found, so the entire observation was 
used. Because the source is 5' off axis, the point spread 
function was fitted with an ellipsoidal Gaussian (1.9" and 
1.1" along the two axes, rms values). Also, the pixel with 
the highest number of counts is offset by 0.8" from the 
center of the fitted ellipse. Taking these uncertainties into 
account, we conservatively estimate a positional error of 
0.7" (1-c). The final Chandra position is: a = 03h 18m 
22.27s±0.12s, S = -66° 36' 03.8"±0.7". 

In order to check the accuracy of the Chandra aspect 
solution, we exploited the presence in the field of view of a 
quite peculiar supernova, SN 1978K, that shows powerful 
radio and X-ray emission. The Chandra position of SN 
1978K is a = 03h 17m 38.69s, S = -66° 33' 03.6" (J2000), 
within 0.46" from the accurate (0.1") radio position of 
Ryder et al. (1993). This is consistent with the expected 
Chandra aspect accuracy. 

2 Based on observations collected at the European Southern Observatory, 



The position of NGC 1313 X-2 was previously deter- 
mined from the ROSAT HRI (Stocke et al. 1995; Schlegel 
et al. 2000) and XMM EPIC-MOS (Miller et al. 2003) im- 
ages. Typical 1-a error boxes are <~ 3" for ROSAT HRI 
and - 2" for XMM EPIC-MOS. The ROSAT and XMM 
positions and corresponding error boxes are summarized 
in Table 1. 

2.2. X-ray spectrum and lightcurve 

We analyzed the XMM EPIC data from both the MOS 
and PN cameras (operated with the medium filter). The 
EPIC-PN spectrum is reported here for the first time and 
was extracted directly from the observation data file be- 
cause the automatic pipe-line processing failed to produce 
a standard event list for the EPIC-PN camera. Both re- 
duction procedures (epchain and epproc) were used to 
extract the data obtaining similar results (differing typi- 
cally by a few percents). Data screening, region selection 
and event extraction were performed using standard soft- 
ware (XMMSELECT v 2.43.2). An analysis of the MOS 
and PN light curves shows that solar flares are present 
in both datasets. They were filtered out using the stan- 
dard criterion (total off-source count rate above 10 kcV < 
5 counts s _1 for MOS and < 15 counts s _1 for PN). We 
extracted the source counts from a circle of 40" and 30" 
for the MOS and PN cameras, respectively. The proxim- 
ity of the source to one of the CCD edges in the EPIC- 
MOS data requires some care. We eliminated the area of 
a box aligned and superimposed to the CCD boundary to 
avoid contamination from bad pixels close to the source. 
The background was selected from a circle of 60" in a 
nearby source- free region of the same CCD. Ancillary and 
response files were produced using the appropriate XMM- 
SELECT tasks. Data were grouped to require at least 20 
counts per bin for the MOS data and 40 counts per bin for 
the PN data, and were then analyzed and compared with 
different models using XSPEC (v 11.2.0). To minimize the 
effects of possible relative calibration uncertainties, the fit 
of the MOS1, MOS2 and PN spectra were performed with 
an overall normalization constant (those of the two MOS 
cameras differ by ~ 10%, while that of the PN instrument 
is larger by ~ 25%). The count rates are 0.08 counts s _1 
for the MOS cameras and 0.25 counts s _1 for the PN. 

In order to reconstruct the X-ray variability history of 
NGC 1313 X-2, we have carefully re-analyzed also the 
ROSAT and ASCA observations. Extraction regions for 
the ASCA SIS data were chosen with care to avoid con- 
tamination from the CCD edges and SN 1978K. Spectra 
were grouped to require at least 15 counts per bin for the 
ROSAT data, and 20 counts per bin for the ASCA data. 

The results of the spectral analysis are listed in Table 2. 
The statistics of the XMM EPIC spectrum is significantly 
improved including the EPIC-PN data (the PN camera 
has almost twice more counts than each single MOS in- 
strument). An absorbed power-law does not provide a 
satisfactory fit of the joint MOS1, MOS2 and PN spec- 
tra (xled = l-^ f° r 249 d.o.f.). Two components models 
provide a significant improvement over single component 
ones. The best fit is obtained with an absorbed soft, ther- 
mal component plus a power-law. Adding a MCD model to 
the power-law results in an improvement of the fit which is 

Chile, Program number 68.B-0083(A). 
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significant at the ~ 4.5/5 (epchain/epproc) a level. Fig- 
ure 2 shows the results for a MCD+power-law fit. The 
resulting best fitting parameters are kT = 200±^° eV, 

r = 2.231^09 and n h = 3.13j£j? x 10 21 cm- 2 for the 
inner disk temperature, photon index and column density, 
respectively (see Table 2). There are residuals in the fit 
(mainly in the EPIC-PN spectrum) that suggest the pos- 
sible presence of emission lines. We emphasize that the 
EPIC-PN data provide marginal evidence for the presence 
of a soft component even at low metallicities. Reducing 
the abundances of the absorbing gas at 0.5 solar, a simple 
power-law fit of the EPIC-PN data has xL d = 1-24 (89 
d.o.f), while a MCD+power-law fit gives Xred = 1-15 (87 
d.o.f.). Finally, it is worth noting that the value of kT is 
~ 25% larger (although consistent within l-er) than that 
derived by Miller et al. (2003). 

In Figure 1 we report the X-ray flux derived from all the 
available observations of NGC1313 X-2. The fluxes were 
consistently derived from the best fit parameters of the 
X-ray spectral analysis reported in Table 2. An approxi- 
mate estimate of the errors, based on counting statistics, 
is 5-10%. The Chandra point is not included in Figure 
1 because of pile-up problems. The unabsorbed 0.2-10 
keV flux from both the XMM EPIC-MOS and PN instru- 
ments agree within ~ 20% and give an average value of 
2.4 x 10~ 12 erg cm _2 s _1 . This value is lower by a factor 
~ 2 than that estimated by Miller et al. (2003). Variabil- 
ity of up to a factor 2 on a timescale of months is clearly 
present and it is reminiscent of the behavior observed in 
Galactic X-ray binaries. If uncertainties in the best-fitting 
spectral parameters are taken into account, the amplitude 
of variability is reduced but not eliminated. This suggests 
that a compact object is present in NGC 1313 X-2. If the 
emission is isotropic and the distance is that of the host 
galaxy (~ 3.7 Mpc; Tully 1988), the X-ray luminosity in 
the 0.2-10 keV range is L x ^ (3-6± 0.5) x 10 39 erg s" 1 . 
If at maximum the source radiates at the Eddington limit 
LEdd, the BH mass is ~ 50 M@. Sub-Eddington accretion 
would imply an even larger mass. 

Interestingly, comparing data from the first and second 
epoch ASCA observations, the X-ray flux of NGC 1313 
X-2 appears to increase with increasing spectral hardness 
(see Table 2). This behavior is similar to that observed in 
the ULXs of the Antennae galaxy (Fabbiano et al. 2003) 
and is opposite to what is usually seen in Galactic BH 
X-ray binaries (e.g. Cyg X-l). 

3. OPTICAL OBSERVATIONS 

Optical images of the field of NGC 1313 X-2 in the R- 
band (Bessel filter) were taken on 16 January 2002 with 
the 3.6 m telescope of the European Southern Observa- 
tory (ESO) at La Silla (Chile). We used EFOSC2 with a 
Loral/Lesser CCD of 2048x2048 pixels yielding a field of 
view of ~ 5' x 5' at a resolution of 0.314"/pixel (re-binned 
by a factor 2). The night was clear with a seeing of about 
1". Four images were obtained for a total exposure time 
of 1320 s. Standard reduction of the data (including bias 
subtraction and flat-field correction) was performed within 
the IRAF (v 2.12) environment. 

A spectrum of one of the field objects (object A; see 
below) was secured on the same night. We performed low- 
resolution (13.4 A, grism#4) spectroscopy for a total expo- 



sure time of 1200 s. After applying standard corrections 
and sky subtraction, cosmic rays were removed and the 
spectrum was corrected for atmospheric extinction. At the 
time of the optical observations object A had already been 
imaged at the 1.1 m Las Campanas telescope by Stocke et 
al. (1995) and was considered a possible counterpart of 
NGC 1313 X-2. Although the new accurate Chandra posi- 
tion rules out an association with this object (see below), 
the spectrum can be used to gain insight on the proper- 
ties of a surrounding nebula, possibly associated with the 
X-ray source (see e.g. Pakull & Mirioni 2002). 

3.1. Astrometry and photometry of field objects 

Our four ESO images were astrometrically calibrated 
using an IRAF task (PLTSOL) and performing a poly- 
nomial interpolation starting from the positions of GSC2 
ESO field stars. The internal accuracy of this procedure 
was estimated comparing the actual positions of a number 
of GSC2 stars not used for astrometric calibration with 
the positions contained in the catalog. The accuracy is 
0.3" (1-c). The four calibrated images were then summed 
together and the resulting image is shown in Figure 3. 

In order to check for the relative systematics between 
the optical and X-ray astrometric calibrations, we used 
again the position of SN 1978K. This supernova is inside 
the Chandra field of view but outside our optical image. 
Thus, we analyzed also an archival image of SN 1978K 
(from the Padova-Asiago Supernova Archive) taken on 13 
September 1999 with the same telescope and a similar in- 
strumental set-up (ESO 3.6m+EFOSC/2.9+R#642, ex- 
posure time 180 s). After calibrating the archival image, 
the position of SN 1978K is a = 03h 17m 38.605s, S = -66° 
33' 03.13" (J2000). This is within 0.28" from the radio po- 
sition of Ryder et al. (1993), improving significantly upon 
the previous optical position by the same authors. The 
difference between the centroids of the optical and Chan- 
dra positions of SN 1978K is 0.69" (a opt - a x = -0.085s, 
S op t — o~x = -0.47"). Although this difference is small and 
comparable with the statistical errors, we decided to ap- 
ply this correction to the Chandra position of NGC 1313 
X-2 to eliminate any systematic error between the optical 
and X-ray astrometric calibrations. The resulting Chandra 
position of NGC 1313 X-2 is reported in Table 1. 

The photometry of the objects in our optical image was 
performed calibrating the frame with the i?-band mag- 
nitudes of 23 stars from the SuperCosmos Sky Survey 
(Hambly, Irwin & MacGillivray 2001) homogeneously dis- 
tributed over the field of view. The internal accuracy of 
this calibration is 0.2 mag. Aperture (5" radius) magni- 
tudes are reported in Table 1. 

3.2. Spectroscopy of the emission nebula 

The two-dimensional spectrum of the field around ob- 
ject A (Figure 4) shows clear emission lines extending for 
tens of arcsecs from east to west, confirming the existence 
of an extended (<~ 400 pc) optical emission nebula that 
was first found in deep H a images by Pakull & Mirioni 
(2002). A one-dimensional spectrum of the field was ex- 
tracted over an aperture of 0.9" (3 pixels) from two regions 
east-ward and west-ward of the position of object A and 
adjacent to it (Figure 5). Wavelength and relative flux 
calibration were applied to the data. The one-dimensional 
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spectrum shows strong emission lines of H Q , H^, [SII] AA 
6717-6731 A, [01] A 6300 A and [OIII] AA 4959-5007 A. 
The shift of the centroid of the lines (<~ 10 A) is consistent 
with the recession velocity of the galaxy and indicates that 
the emission nebula is located in NGC 1313. 

It is worth emphasizing the abrupt change in the abso- 
lute and relative intensity of the emission lines from east 
to west, indicating variations in the physical conditions 
and/or geometry of the emission nebula. In particular, 
strong emission from [OIII] is present on the east side but 
almost absent on the west side, while emission from H a , 
[SII], [OI] and other elements is present on the west side 
but weaker or absent on the east side. A careful inspection 
of the line intensity profiles (in particular those of H Q and 
[SII]) reveals a fairly symmetric, broadly peaked profile, 
centered at <~ 2" west of the position of object A, and a 
weaker, roughly constant intensity component extending 
in the east direction. 

4. DISCUSSION 

Our Chandra position of NGC 1313 X-2 (Table 1) 
is shown in Figure 3, together with the ROSAT HRI 
(Schlegel et al. 2000) and XMM EPIC-MOS (Miller et 
al. 2003) error boxes, overlaid on our ESO image. All 
measurements are consistent within 1-cr. The distance of 
the centroids of objects A, B and D with respect to the 
Chandra position is 3.6", 4.1" and 7.3", respectively. Even 
taking into account the statistical error on the optical po- 
sitions (0.3"), these three objects can be ruled out at a 
significance level of at least 3-ct. On the other hand, ob- 
ject C is inside the Chandra error box and its position 
coincides within 1-cr with that of NGC 1313 X-2, making 
it a likely counterpart. 

From the maximum absorbed X-ray flux of NGC 1313 
X-2 ( fx ~ 2 x 10~ 12 erg cm~ 2 s -1 ) and optical magnitude 
of object C (R = 21.6), we estimate fx/ f op t ~ 500. This 
value is very high, in agreement with the suggestion by 
Cagnoni et al. (2002) that ULXs can be selected on the 
basis of their large fx/ f opt- Only Isolated Neutron Stars 
(INSs), heavily obscured AGNs and luminous X-ray bina- 
ries can reach such large values of the X-ray/optical flux 
ratio. INSs are extreme in this respect, with B w 25 op- 
tical counterparts and typical X-ray-to-optical flux ratios 
> 10 5 (see e.g. Kaplan, Kulkarni & van Kerkwijk 2003) 
The presence of the relatively bright (R ~ 21.6) object C 
in the Chandra error box makes this possibility unlikely. 
Furthermore, known INSs exhibit different spectral prop- 
erties with no significant variability. On the other hand, a 
heavily obscured AGN is expected to have a rather hard X- 
ray spectrum and to emit significantly in the near-infrared 
(see e.g. Brusa et al. 2002). Given the X-ray luminosity of 
NGC 1313 X-2, an infrared magnitude K w 12 is expected 
were it an obscured AGN. The lack of any IR counterpart 
on a K image of the 2MASS All Sky Image Service down to 
a limiting magnitude K ~ 14 (10-cr) and the softer X-ray 
spectrum of NGC 1313 X-2 make this possibility unlikely, 
although a low resolution optical spectrum of object C is 
definitely required to settle this issue. Our accurate Chan- 
dra position and optical identification favor a very lumi- 
nous X-ray binary in NGC 1313 as the likely explanation 
for NGC 1313 X-2. As a reference, the X-ray/optical flux 
ratio of persistent BH binaries at maximum is > 10 — 100, 



while that of soft X-ray transients in outburst can reach 
2000 (Masetti 1997). This is in line with the alleged bi- 
nary nature of ULXs and is consistent with the observed 
properties of this source, such as the X-ray variability and 
the observed X-ray spectrum, including the presence of a 
soft component probably produced by an accretion disk. 

If indeed NGC 1313 X-2 is a black hole binary, the X- 
ray spectral parameters, in particular the temperature of 
the MCD fit (hereafter referred to as Tmcd), can be used 
to estimate the BH mass. This is similar to what done by 
Miller et al. (2003) using the normalization of the MCD fit 
and, as discussed below, we reach similar conclusions. The 
effective temperature of a standard accretion disk depends 
on radius as T 4 = (3GM BH M /87rerr 3 „)(r m /r) 3 / 4 [l - 
(rin/r) 1 / 2 ], where M is the accretion rate and r in is the in- 
nermost disk radius (e.g. Frank, King & Raine 1980). As- 
suming that Tmc d represents an estimate of the maximum 
disk temperature, it is (SGMbhM /Snarf^ 1 / 4 = uTmcd, 
with a ~ 2. Neglecting relativistic corrections and as- 
suming that the disk terminates at the innermost stable 
circular orbit of a Schwarzschild BH, it is: Mbh/Mq = 
(Mc 2 lL Edd )f\uT MCD l\.h x 10 7 K)- 4 , where / is a color 
correction factor (/ ~ 1.6 — 1.7, Shimura & Takahara 1995; 
Zampieri, Turolla & Szuszkicwicz 2001). Given the strong 
dependence of Mbh on temperature, any uncertainty in 
the accretion physics and radiative transfer may induce 
significant errors in the resulting value of Mbh- So, the 
inferred spectroscopic measurement of the BH mass should 
be taken simply as an approximate estimate. The low in- 
ner disk temperature obtained from the two-components 
fit to the XMM EPIC spectrum (kT ~ 200 cV) implies 
M BH ~ 90/ 4 M (a ~ 2) for Eddington limited accre- 
tion, somewhat larger than that derived from the flux. 
This result removes the need for a rapidly spinning BH 
(invoked by Makishima et al. 2000 from an analysis of the 
ASCA data) and agrees with the conclusion of Miller et al. 
(2003) that NGC 1313 X-2 contains an intermediate mass 
BH. The large inferred BH mass does not require beamed 
emission. Then, the estimated accretion rate (assuming 
10% efficiency) is M ~ 10 _7 M Q yr _1 , forcing the mass 
reservoir to be a companion star. 

From the apparent magnitude (R — 21.6) and absorp- 
tion (Ar ~ 1.7, computed from the X-ray best fitting col- 
umn density N H ~ 3 x 10 21 cur 2 ; Bohfin et al. 1978) of 
object C, we estimate an absolute magnitude Mr ~ —7.9 
and a luminosity in the range ~ 7 x 10 — 1O 6 L0, depend- 
ing on the adopted bolometric correction. If this originates 
from the companion star, the inferred luminosity is consis- 
tent with a « 20M Q main sequence star ora~ 15 — 20M Q 
evolved OB supcrgiant (e.g. Bowers & Deeming 1984), 
making NGC 1313 X-2 a high-mass X-ray binary. In lu- 
minous Galactic X-ray binaries, the reprocessed optical 
emission from the disk may be significant. Assuming that 
20-30% of the X-ray flux produced in the innermost part 
of the accretion disk intercepts the outer regions, for real- 
istic values of the albedo (> 0.9; e.g. de Jong, van Paradijs 
& Augustcijn 1996) few percents of the X-ray luminosity 
(« 10 38 erg s _1 ) can be absorbed and re-emitted in the op- 
tical band. Characteristic emission lines of X-ray-ionized 
H, He or N, typically seen in luminous Galactic X-ray bi- 
naries should then be detectable in the optical spectrum. 
Also X-ray heating of the companion star itself may con- 
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tribute to the optical emission. If the optical luminosity 
comes in part from X-ray re-irradiation, the mass of the 
companion would be lower. Taking M 2 ~ 2OM as an 
upper limit for the mass of the companion, the mass ratio 
of the binary is q = M 2 /M BH < 0.4/~ 4 < 1. Writ- 
ing the binary separation as a — 2.16i?2[<z/(l + q)]^ 1 ^ 3 
(Paczyhski 1971), the orbital period of the system is 
P > 0.15(R 2 /100RQ) 3/2 f 2 (M BH /50M & )~ 1 / 2 yr. Ac- 
cording to King et al. (2001), the system should not be 
a persistent X-ray source. 

Although difficult to reconcile with the properties of the 
environment surrounding NGC 1313 X-2 (see below), we 
can not rule out also that the optical emission detected 
in the Chandra error box originates from a stellar cluster 
(see e.g. the case of a ULX in NGC 4565; Wu et al. 2002), 
in which case NGC 1313 X-2 may be a low-mass X-ray 
binary in the cluster. 

The mass accretion rate required to produce the ob- 
served luminosity may in principle be provided by Roche- 
lobe overflow from an evolved companion or by a wind 
from a supergiant. In the first case, evolutionary swelling 
of the companion keeps pace with the increase in Roche 
lobe size and the system remains self-sustained: accretion 
is likely to proceed through a disk. In the second case, 
assuming 10% accretion efficiency and that the BH can 
capture ~ 1% of the mass outflow, the wind must be very 
powerful (M ~ 10 _5 M Q yr _1 ). A lower efficiency would 
require too high a gas supply, so a disk is needed even 
in a wind-fed system. In this case, however, the disk is 
probably much smaller than in a Roche-lobe overflow sys- 
tem and the optical emission dominated by the supergiant. 
On the other hand, in a Roche-lobe overflow system, an 
extended, possibly re-irradiated accretion disk should con- 
tribute heavily in the UV and B bands, producing strong 
emission lines. Thus, the two modes of mass supply are 
likely to be distinguishable by optical spectroscopy. 

We now turn to discuss how our optical observations 
can be used to constrain the environment of NGC 1313 
X-2. Figures 4 and 5 reveal that NGC 1313 X-2 is likely 
to be associated with an optical emission nebula, recogniz- 
able also in a H Q image taken by Pakull & Mirioni (2002). 
From the velocity (80 km s -1 ) and flux of fig, they derive 
an impressive mechanical energy of 3 — 10 x 10 52 erg for 
the expanding ionized gas, and suggest that the nebula 
is inflated by a relativistic jet from NGC 1313 X-2. Our 
measured ratio of [SII]/H Q (~ 0.5) is consistent with that 
expected from a shock- ionized supernova remnant, a stel- 
lar wind-shocked nebula or diffuse ionized gas (Matonick 
& Fesen 1997). However, the inferred diameter and energy 
of the nebula are too large to be consistent with a single su- 
pernova event, unless it was produced by a hypernova sim- 
ilar to SN 1998bw (see e.g. Iwamoto et al. 1998). In fact, 
it could be the result of several explosion events (multiple 
supernova remnant) or be originated by the intense wind 
of hot stars, possibly the parent stellar association of NGC 
1313 X-2. As discussed in the previous section, the nebula 
appears to have some internal structure: a comparatively 
brighter, fairly symmetric component west of the position 
of NGC 1313 X-2 and a weaker, slightly elongated one ex- 
tending in the east direction. The brighter part of the neb- 
ula has [SII]/H a =0.58, the weaker one has [SII]/H Q =0.44 
and intense [OIII] emission. Different possibilities may ex- 



plain the irregular appearance of the nebula. As suggested 
by Pakull & Mirioni (2002), the varying line intensity may 
be caused by reprocessed emission from the X-ray ionized 
interstellar medium where the physical conditions (in par- 
ticular the density) vary on a scale ~ 100 pc. However, 
the nebular emission may also arise from two physically 
distict components: a wind-shocked nebula produced by a 
possible parent stellar association of NGC 1313 X-2 and 
a multiple supernova remnant. This hypothesis seems to 
be confirmed also by the marginal detection of (possibly 
extended) UV emission in an image of the XMM Optical 
Monitor (see Figure 6), in coincidence with the brighter 
component. Clearly, the weaker component may still be 
a jet-inflated nebula, as suggested by Pakull & Mirioni 
(2002). Finally, we note that, although NGC 1313 X-2 is 
somewhat hotter and much more luminous, the [OIII] sig- 
nature in the eastern portion of the nebula is reminiscent 
of predictions for the radiation-limited nebulae around su- 
pcrsoft sources (Di Stefano, Paerels & Rappaport 1995; 
Chiang & Rappaport 1996). 

It is interesting to note that object A, the possible as- 
sociation of which with NGC 1313 X-2 was discussed by 
Stocke et al. (1995), lies very close both to the Chan- 
dra error box and the point where the intensity of the 
nebular emission lines suddenly changes. The continuum 
spectrum of this object, obtained after subtracting off the 
emission line spectrum of the nebula in the adjacent re- 
gions, was compared to template stellar spectra (Pickles 
1998) and turns out to be in fair quantitative agreement 
with that of a G-M supergiant, but not with that of late- 
type dwarfs. This result is independent of reddening. The 
absolute magnitude of object A is Mr ~ —9.7, and the 
luminosity L ~ 5 x 10 5 L Q (assuming a bolometric correc- 
tion appropriate for a K star). Object A may be a very 
massive (~ 30M Q ) G-M supergiant of radius ~ 1000i? Q 
or a cluster in NGC 1313. The first interpretation would 
support the conclusion that the region in which NGC 1313 
X-2 is located is an active star forming region, in which the 
initial mass function is top-heavy. On the other hand, the 
second possibility appears more likely because ~ 30M Q 
massive stars are extremely rare. However, the red color 
would indicate a rather evolved cluster that would be pro- 
jected by chance on the active star forming environment 
in which NGC 1313 X-2 appears to be embedded. 

A crucial question is how a binary system containing an 
intermediate mass BH may have formed (see e.g. van der 
Marel 2003). The BH progenitor must have been rather 
massive. This is consistent with the fact that NGC 1313 
is likely to have lower than solar metallicity (Z ~ 0.5; 
Zaritsky, Kennicutt & Huchra 1994) and hence mass loss 
was less intense. Such a massive BH may have formed 
through direct collapse without producing a supernova. In 
this way, if the system was born as a binary, it may have 
survived after the collapse of the primary. Although less 
likely, it is also possible that the companion might have 
been captured from a nearby stellar association. In this 
case, it is not possible to exclude that the BH may have 
formed from an early episode of star formation (population 
III). 

It is worth noting that, although the large BH mass does 
not require that the emission is beamed, we cannot rule 
out that a moderate jet activity, producing radio emission 
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(and possibly inflating the emission nebula), is present in 
NGC 1313 X-2 (see e.g. the case of an ULX in NGC 5408; 
Kaaret et al. 2003). However, presently available radio 
images of the field of NGC 1313 X-2 (Sydney University 
Molonglo Sky Survey at 843 MHz and Australia Telescope 
Array at ~ 5 GHz; Stocke et al. 1995) are not sufficiently 
deep to allow detection. 

Optical spectroscopy of object C, narrow band imaging 
of the nebula, deep radio observations and an analysis of 
the short timescale X-ray variability are essential to bet- 
ter assess the physical properties of NGC 1313 X-2 and 
its environment. In particular, even a moderate resolution 
spectrum of object C will make it possible to detect any 
characteristic absorption and/or emission line, and then 
determine its properties and redshift. These observations 
will allow us to strengthen the identification of NGC 1313 
X-2 with an intermediate mass BH and foster our under- 



standing of ULXs. 
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Table 1 

Positions of NGC 1313 X-2 and positions and optical magnitudes of field objects 



Observatory/Instr. Object RA[J2000] DEC[J2000] R magnitude (Bessel-Cousins) Ref. 



ROSATpSRl 


NGC 


1313 


X-2 


03 


18 


22.00±0.50 


-66 


36 


02.3±3.0 




Schlegel et al. (2000) 


XMM/EPIC-MOS 


NGC 


1313 


X-2 


03 


18 


22.34±0.33 


-66 


36 


03.7±2.0 




Miller ct al. (2003) 


Chandra/ ACIS-S 


NGC 


1313 


X-2 


03 


18 


22.18±0.12 


-66 


36 


03.3±0.7 




this work 


ESO/3.6m 


A 






03 


18 


21.97±0.05 


-66 


36 


06.5±0.3 


19.8±0.2 


this work 


ESO/3.6m 


B 






03 


18 


21.56±0.05 


-66 


36 


00.9±0.3 


20.7±0.2 


this work 


ESO/3.6m 


C 






03 


18 


22.34±0.05 


-66 


36 


03.7±0.3 


21.6±0.2 


this work 


ESO/3.6m 


D 






03 


18 


20.96±0.05 


-66 


36 


03.7±0.3 


17.8±0.2 


this work 



a See Figure 3. 



Table 2 

Parameters of the fit of ROSAT, ASCA and XMM observations 



Observatory /Instr. 


Obs. Id. 


Model 


(10 21 cm" 2 ) 


Parameters 


x'ed (dof) 


Fx 

(10~ 12 erg cm" 2 s" 1 ) 


r 0.2-10 koV 

(10~ 12 erg cm" 2 s" 1 ) 


ROSAT/PSPC 


rp600045n00 


Power-law 


1 06+ 214 

± ' UD — 0.50 


r - 2 40+ 1 - 40 


0.83 (20) 


1.4 


2.6 






Blackbody 


17+ - 20 

u - 1 ' -0.15 


kT = 0.27±S;g| 


1.01 (20) 


0.34 


0.38 


ROSAT/PSPC 


rp600504n00 


Power-law 


I 7 4+2-95 
L - '^-1.31 


r = 2.20±J;|° 


1.43 (22) 


1.5 


1.8 






Blackbody 


09 +0 " 23 


kT = 0.35+°°= 


1.32 (22) 


0.37 


0.44 


ASCA/GIS+SIS 


60028000 


Power-law 


4 24+ - 41 


r - 1 99+ 03 


0.99 (171) 


1.3 


3.2 






MCD 


i i r + 0.35 
1 ' lo -0.34 


kT = 1.38+°° 7 7 


0.89 (171) 


2.1 


2.6 


ASCA/GIS+SIS 


93010000 


Power-law 


4 74+0-50 

*■ '*-0.50 


r - 2 43+ 09 


1.28 (139) 


0.50 


2.2 






MCD 


1 06+ - 31 

1 - UD — 0.31 


kT = 1.08+. , , 4 


1.31 (139) 


0.86 


1.0 


XMM/EPIC-MOS+PN 


0106860101 


Power-law 


2 69 +0 18 


r - 2 42+° 07 


1.20 (249) 


2.1 








MCD 


70+ 06 

u - ' u -0.06 


kT = 0.90±° ° 2 2 


2.29 (249) 


0.86 








CompTT c 


1 67+ ' 43 

iu ' -0.45 


kT = S.OOtiH 


1.11 (248) 


1.2 










Jo — U.10_„o3 

r = 4.69lS;iS 












MCD+Power-law 


o -,0 + 0.92 
°- lo -0-37 


fcT = o. 20 j:° : 4 


1.09 (247) 


2.4 










p n Qo + 0.15 

i — z.zo_ 09 












MCD+Power-law d 


o C7+0-94 
a0 1 -0.52 


fcT = 0.20i°' 07 


1.07 (247) 


2.0 










r = 2.23l^; 2 









a Unabsorbed flux in the 0.1-2.0 kcV (ROSAT), 2.0-10.0 kcV (ASCA) and 0.2-10.0 keV (XMM) energy bands 
b Unabsorbed flux extrapolated in the 0.2-10.0 keV band using the web interface to PIMMS (v 3.3) 
c Thcrmal comptonization model with Wicn soft photon input 
d Abundancc 0.5 solar 
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Fig. 1. — The 0.2-10 keV light curve of NGC 1313 X-2 from all presently available data. Fluxes were de-absorbed and, when necessary, 
extrapolated in the 0.2-10 keV interval using the web interface to PIMMS (see Table 2). For the Einstein IPC observation, we adopt the 
value reported by Gioia et al. (1990). For the ROSAT PSPC and ASCA GIS+SIS observations, the fluxes were derived from the best fit 
parameters of the power-law model reported in Table 2. For the ROSAT HRI data (Schlegel et al. 2000) a power-law spectrum with T = 2 
and Nj] = 3 X 10 21 cm -2 was assumed (in agreement with the best fit parameters derived from the spectral analysis). For the XMM EPIC 
data, the flux is calculated from the best fit MCD+power-law model reported in Table 2. 



Ultraluminous X-ray Sources: NGC 1313 X-2 
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[A 9 >1 s ^uuo/suo^oqj] 

Fig. 2. — X-ray spectrum of NGC 1313 X-2 from the XMM EPIC-MOS and PN cameras. The solid line represents the combined best 
fitting model spectrum, while the dashed and dotted lines are the MCD and power-law components respectively. 
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Fig. 3.— ESO 3.6m fl-band (Bessel filter) image of the field of NGC 1313 X-2. The circles show the ROSAT HRI, XMM EPIC-MOS and 
Chandra ACIS-S positions. The estimated 90% confidence radii are 6" for HRI, 4" for EPIC-MOS and 1.4" for ACIS-S. Labels A, B, C and 
D mark the four field objects inside or close to the X-ray error boxes. The insert at the bottom-left shows a larger portion of the image with 
the position of the X-ray source (cross). 




Fig. 4. — Two-dimensional spectrum (ESO 3.6m+EFOSC2+grism#4) of the field around object A. The slit (1.2") is oriented in the 
east-west direction. The wavelength intervals are 4500-5300 A (left panel) and 6150—6900 A (right panel). 



Ultraluminous X-ray Sources: NGC 1313 X-2 
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Fig. 5. — One-dimensional spectrum (Fx) of the nebula around NGC 1313 X-2. The lower (upper) spectrum is extracted in a region 
east- ward (west- ward) of the position of object A (see text for details). The symbol <g) marks a residual contamination from an emission line 
of the sky. 
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Fig. 6.— UV-band (UVW1 filter, 1800-3200 A) exposure of the field of NGC 1313 X-2, obtained with the XMM Optical Monitor. The 
circle shows the Chandra position (90% confidence level). A region of possible diffuse emission is visible NW of the Chandra error box. 



